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Abstract. We detect 1.2 mm continuum emission from 
dust in the gravitationally lensed Lyman break galaxy 
MS 1512+36-cB58. Our detected flux is surprisingly low: 
relative to local star burst galaxies, cB58 appears to pro- 
duce somewhat less far-IR emission than its UV redden- 
ing predicts. After comparing several different estimates 
of the source's dust content, we conclude that the appar- 
ent discrepancy is most likely related to uncertainty in its 
UV spectral slope. Alternate scenarios to account for a 
far-IR "deficit" which rely on a high dust temperature or 
differential magnification are less satisfactory. Our result 
underscores one of the risks inherent in characterizing the 
cosmic star formation history from rest-UV data alone. 

Key words: Galaxies: evolution, fundamental param- 
eters, photometry, starburst; Cosmology: observations; 
Submillimeter 



1. Introduction 

Large samples of actively star-forming galaxies at z > 2 
can now be identified through color selection techniques 
which exploit the Lyman break (e.g., Steidel et al. 1999| ). 
In many respects, these sources resemble UV-bright star- 
burst galaxies at z — 0: the two populations have similar 



bolometric surface brightnesses (Meurer et al. 1997), UV- 
through-optical spectral energy distributions (SEDs: e.g., 
Ellingson et al. 1996. Sawicki & Yee 1998), and blueshifted 
inters tellar absorption from outflowing gas (e.g., Pcttini 
et al. 1998 ). The analogy between Lyman break galaxies 
(LBGs) and local starbursts is not exact, however: LBGs 
appear to be forming stars at rates which are "scaled 
up" by their larger physical sizes (Meurer et al. 1997 ). 
At longer wavelengths, moreover, the analogy can barely 
be evaluated, since few direct far-IR and submillimeter 
observations of LBGs exist. Two recent at temp ts have 
been made to fill this gap. Chapman et al. ( 2000 ) detect 
only one of sixteen LBGs at 450 (im and 850 /jm with the 
Submillimeter Common-User Bolometer Array (SCUBA). 
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Peacock et al. (2000) compare HST and SCUBA observa- 



tions of the Hubble Deep Field and achieve a statistical 
detection at 850 /im of the HST sources with photomet- 
ric redshifts, but can make no statement about specific 
objects with spectroscopic redshifts. The limitations of 
these studies motivated us to observe the gravitationally 
len sed s ystem MS 1512+36-cB58 (hereafter cB58: Yee et 
al. 1996). Thanks to magnification by the z = 0.37 cluster 
MS 1512+36, this z = 2.7 LBG was already well-studied 



in the rest UV and optical (e.g., Ellingson et al. 1996; Pct- 
tini et al. 2000. Teplitz et al. 2000), and promised to be 



unusually detectable at longer wavelengths. 
2. Observations 

From 3-7 March 2000, we observed cB58 in on-off photo- 
metric mode with the 37-element Max-Plan ck Mi llimeter 
Bolometer (MAMBO) array (Kreysa et al. |l998j ) at the 



IRAM 30m telescope. We pointed the array's central pixel 
at the optical position (J2000 coordinates 15:1 4:22.2 and 
+36:36:24.4) determined by Abraham et al. ( |L998[ ). We 
chopped the secondary mirror by 53" in azimuth at 0.5 Hz 
and nodded by the same 53" throw every ten seconds; the 
HPBW of our beam was 11" at 1.2 mm. We reduced the 
data with the MOPSI package, using all pixels in the first 
ring around the central pixel for estimating and removing 
correlated fluctuations due to the sky background. Our 
flux scale of 6250 counts Jy _1 was based on observations 
of planets. If we include all data from our total on+off in- 
tegration time of 200 minutes, we achieve a 6. Oct detection 
of 1.32+0.22 mJy. For the purposes of this paper, however, 
we will exclude 60 minutes of data from our last observing 
night, for which pointing and sky background instabilities 
coupled with a large opacity correction (r ~ 0.3 vs. < 0.16 
over the rest of our run) reduce our confidence in the pho- 
tometry. Including a 10% uncertainty in our flux scale then 
yields a final [4.4ct] detection of 1.06 ± 0.35 mJy. 

3. A comparison with local starbursts 



Heckman et al. ( |1998| ) and Meurer et al. ( |1999[ ) have 
showed that local starburst and star-forming galaxies obey 
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scaling relations between various global properties and the 
parameters of their UV spectra. We can therefore test 
the proposition that cB58 is a "scaled-up" local starburst 
by seeing if it fits those relations which are independent 
of distance (and lensing magnification). A key parameter 
here is the UV spectral slope (3 (fx cx A' 3 ), whose ob- 
served value indicates the amount of reddening by dust 
(i.e., from some intrinsically blue /3rj) for an assumed star 
formation history. Because the same dust which reddens 
stellar UV continua will re-radiate the energy which it ab- 
sorbs, it is natural that UV-bright local systems should 
obey a relationship between (3 an d the ratio of far-IR to 
UV emission. Meurer et al. ( 1999 ) determine this relation 
for a sample which includes only galaxies whose starbursts 
are small enough to fit entirely within the IUE aperture. 
They quantify the fraction of reprocessed UV emission as 
a ratio between 

F Fm = 1.26xl(r 11 [2.58/„(60) + / I ,(100)] ergcm" 2 s -1 (l) 

for 60 /xm and 100 /Ltm IRAS flux densities in Jy, and 
-P1600 = v fy evaluated at v = c/1600A. Figure |l| plots 
this ratio versus (3 for the 43 sources which (a) appear in 
their Figure 1, and (b) have flux densities in the Version 
2.0 IRAS Faint Source and/or Point Source Catalogs. The 
local starbursts occupy a predictable locus in the diagram, 
extending from blue, metal-poor systems at the lower left 
to reddened, dusty systems at the upper right. 

To place cB58 on this diagram, we first measu re (3 
from the rest-UV spectrum of Pettini et al. ( |2000 ), ac- 
cording to the prescription of Meurer et al. ( 1999|) for 
the sake of consistency. A least-squares fit (with iterative 
±3tr rejection) to eight spectral windows between 1100 
and 1975 A gives (3 = -1.12 ± 0.10, slightly bluer than 
the -0. 8 (f3 = —2 — a for /„ cx v a ) measured by Pettini 
et al. ( 2000 ) from their fit to co ntinuu m flux densities at 
1300 and 1800 A. Meurer et al. fll999| ) also add an offset 
A/3 = —0.16 ± 0.04 for each source lacking spectral cover- 
age longwards of 1975 A; making the same adjustment for 
cB58 yields our final value of = -1.28 ± 0.14. Next, we 
can estimate Fi 6 oo f rom t he Vab = 20.64 ± 0.12 observed 
by Ellingson et al. (1996), corresponding to a rest wave- 



length of 1474.9 A (for z = 2.7290: Teplitz et al. |2000|) 

7 i474.9 = vfv (observed 5500 A) 

= 1.10 ± 0.13 x 10~ 13 erg cm~ 



2 s" 1 



Scaling this flux by a factor (1600/1474.9) 1+/3 = 0.98 ± 
0.01 then gives Fi 600 = 1.07 ± 0.14 x 10~ 13 ergcm~ 2 s" 1 . 

Our final step is to derive Ffir from our observed 
1.2 mm flux density. We must multiply the observed uf v — 
2.65 x 10~ 15 erg cm -2 s _1 by a rest-frame ratio of the form 



J FIR. 
Vfv 



1.26 x 10~ n 



2.58 



/„(60) /„(100) 



/„(321.8) /„(321.8) 



(2) 



where now v = c(l + z)/1.2mm = c/321.8/im. We must 
evaluate the term in brackets by assuming some rest-frame 
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Fig. 1. Far-IR-to-UV flux ratio vs. UV spectral slope. Es- 
timates for local starbursts are triangles (for IRAS PSC 
flux densities), circles (for IRAS FSC flux densities), or 
pairs of arrows (for FSC 60 /im detections and 100 /im up- 
per limits) . The fit and (for NGC 2537) representative er- 
ror bars are taken from Meurer et al. (1999). cB58 is an 
open square. 



SED for cB58. Since comparably deep SCUBA observa- 
tions of cB58 at 850 /im have yielde d both a detection of 
4.2 ± 0.9 mJy (van der Werf etal. |200l[ ) and a 3ct up- 
per limit of 3.9 mJy (Sawicki [2001 ), the prospects for an 
observationally-derived SED are not encouraging. Instead, 
we simply extend the scaled-up starburst hypothesis for 
LBGs: we assign to cB58 a dust emissivity index e = 2 and 
dust temperature Td = 33.0 ± 4.5 K (the mean and dis- 
persion of the 60/100 /zm color temperatures for the local 
galaxies in Figure |l|, as determined from e = 2 modified 
Planck function fits). We then predict an observed flux 
in the rest-FIR window of 1.28±J;|| x 10~ 13 erg cm" 2 s _1 , 



with large and asymmetric error bars dominated by the 
uncertainty in Td. We emphasize that our extrapolated 
FIR flux is conservatively large: adoption of e = 1 and the 
corresponding mean Td = 40. 7 ±7.1 K would have lowered 
our estimate by a factor of two. We have also ignored the 
potential impact of optically thick dust emission, or of a 
dust component colder than 33.0 K which could dominate 
the 321.8 /mi emission but contribute little to Tfir. 

When plotted in Figure |l], cB58 falls below and/or to 
the right of all objects except UGC 6456. While the scatter 
in the local sample is large, it is s trictly true that relative 
to the fit by Meurer et al. ( 1999| ), cB58 appears to have 
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too red a UV spectral slope for its -Ffir/-Fi600, and too 
low a far-IR-to-UV ratio (by a factor of ~ 4.7) for its f3. 

4. Four estimates of the dust content in cB58 

In physical terms, Figure |l| suggests that cB58 has some- 
what more dust than its far-IR emission predicts or some- 
what less dust than its UV reddening predicts. We can 
explore these two possibilities by making four distinct es- 
timates of the source's dust content on the basis of the 
local starburst analogy. We rely in particular on the effec- 
tive attenuation curve fc'(A) and other empirical relations 
presented by Calzetti ( 2000| ) and references therein. 

I. The far-IR-to-UV ratio 



From their local starburst sample, Meurer et al. (1999) 
derive an empirical relation between A(1600A)~ the ex- 
tinction at 1600 A- and the -Fbtr/ fi600 ratio: 



log(iWFi 60 o) =log(l0°- 



4A(1600 A) 



0.076 



(3) 



From our detection of cB58, we deduce ,4(1600 A) 



0.75 



+0.63. 
-0.47' 



we then use the Calzetti attenuation curve 



to evaluate the color excess of the stellar continuum, 
E S (B -V) = A(X)/k'(X) = 0.075™. 

II. The Ha-to-H/3 ratio 

Because ionizing stars are generally younger and more 
embedded than those which dominate the UV continuum, 
the color excess measured from observations of (gas) re- 
combination lines, E g (B — V), tends to exceed E S (B — V) 
and can be straightforwardly measured using any optical 
extinction curve k(X). In cB58, Teplitz et al. ( |2000|) re- 
port a Balmer decrement of 3. 23 ± .4, which for the LMC 
extinction curve of Howarth ( 1986 ) implies E g (B — V) = 
0.12 ±0.1 2 (rat her than the former authors' — 0.27). Since 
Calzetti (2000) finds that local starbursts obey the pro- 
portionality E S (B - V) = 0.44 E g (B -V), we conclude 
E S (B-V)^ 0.053™. 

III. The UV line spectrum and UV-to-Ha ratio 

The r est-UV spectrum of cB58 obtained by Pettini et 
al. ( 2000 ) is of a quality which permits deta iled popula- 
tion synthesis modelling. Pettini et al. (2000) themselves 
argue that the presence of P Cygni features due to mas- 
sive stars supports a history of continuous star formation, 
with an Initial Mass Function (IMF) which is Salpeter up 
to M ~ 100 M Q . De Mello et al. ( |2000D use the strengths 
of both O and B star features to argue that continuous 
star formation has proceeded for 25-100 Myr, an age range 
whose lower end agrees with the SED fit of Ellingson et 
al. ( |1996| ). For a metallicity Z = 0.4 Z© (comparable to 



the Z ~ 0.32 Z Q measured by Teplitz et al. |2000), the 
appropriate Starburst99 models of Leitherer et al. (1999) 
predict that the intrinsic ratio Rq = -F1475/-FHQ will rise 
from ~ 98 to ~ 120 over the course of the 25-100 Myr in- 
terval. cB58 has an observed R = 90; the fact that R < Ro 
can be attributed to the fact that extinction is higher in 



the UV than in the optical. Using the Calzetti fc'(A), the 
LMC fc(A), and the relation between color excesses, we 
can derive a single equation in E S (B — V): 



log (i?/i? ) = 0AE s (B - V) x 

(fc(6563 A)/0.44 - fc'(1475 A) 



(4) 



whose solution is E S (B - V) = 0.042 ± 0.021. Such low 
values are required to avoid discordant UV and Ha-based 
estimates of the star formation rate in cB58. 

IV. The UV spectral slope 

According to the empirical relation derived by Meurer 
et al. ( |1999D , A(1600 A) = 1.99(/3 - /3 ) for f3 Q = -2.23. 
Using the Calzetti fc'(A), we then derive E S (B — V) = 
0.190 ±0.028. From the analogous relation between (3 and 
E S (B — V) in Calzetti ( [2000 ) and her suggestion of /3 = 
-2.3, we would obtain E S (B - V) = 0.236 ± 0.032. 



5. Conclusions 

Although our calculations of E S (B — V) in cB58 take un- 
certain measurements and apply scaling relations about 
which there is substantial local scatter, it remains strik- 
ing that three of the four estimates arc in excellent agree- 
ment with each other. Only the UV spectral slope gives 
a value which is too high by a factor of 2-3. In principle, 
use of a steeper (e.g., SMC) extinction curve rather than 
the Calzetti fe'(A) would lower the color excess estimated 
from the observed /? — /3o- However, a steeper extinction 
curve would also reduce the values of E s (B — V) derived 
from -Ffir/-Fi600 arid population synthesis models; the dif- 
ference is real. 

The preponderance of evidence for a low dust content 
in cB58 argues against analyses which would describe its 
location in Fi gure |l| solely in term s of a far-IR deficit (e.g., 
Sawicki |200l| ; van der Werf et al. |200l| ). We can erase an 
apparent shortfall in -Ffir/ ^1600, for example, if we pre- 
dict Ffir from vf v (321.8 /im) under the assumption that 
the dust in cB58 is as hot as the dust in the most ex- 
treme of the local starbursts, Tololo 1924-416 (T d = 50.3 K 
for e = 2). However, since none of the other 42 systems 
has Td > 40 K, we would essentially be discarding the lo- 
cal starburst analogy for LBGs while failing to explain 
how two independent lines of evidence (the Ha/H/3 and 
UV/Ha ratios) could still favor a low E S (B — V). Simi- 
larly, we might imagine that the observed far-IR-to-UV 
ratio is depressed from its intrinsic value because UV 
emission reddened by dust along a favored line of sight 
is more highly magnified than the bulk of the far-IR emis- 
sion. This scenario too seems implausible. Since the best 



lensing model of Seitz et al. (1998) predicts that a large 
fraction (~ 0.54 — 0.67) of the background source is mag- 
nified into the cB58 arc, and since UV continuum and Ha 
emission have broadly simi lar sp atial distributions in local 
starbursts (Conselice et al. 2000| ), we would require excep- 
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tional tuning of the lensing geometry to produce agree- 
ment in all estimates of E S (B — V) except that from (3. 

The most straightforward conclusion about cB58 is 
that its offset from the local relation in Figure [I] has more 
to do with (3 than with Ffir/ -Fi600 _ i-e., that the local 
starburst analogy holds more robustly for all other proper- 
ties of cB58 than it does for the UV spectral slope. Obser- 
vational systematics may be partly responsible: a spectro- 
scopic (3 can be artificially reddened by differential refrac- 
tion, while photometry can mislead due to intergalactic 
absorption- as in the case, perhaps, of t he g b and mag- 
nitude for cB58 used by Ellingson et al. ( 1996 ) to derive 



made use of the NASA/IPAC Extragalactic Database (NED), 
which is operated by the Jet Propulsion Laboratory, Califor- 
nia Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 
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